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ABSTRACT 

The South Pole Telescope (SPT) is currently surveying 2500 deg^ of the southern sky to detect 
massive galaxy clusters out to the epoch of their formation using the Sunyaev-Zel'dovich (SZ) effect. 
This paper presents a catalog of the 26 most significant SZ cluster detections in the full survey region. 
The catalog includes 14 clusters which have been previously identified and 12 that are new discoveries. 
These clusters were identified in fields observed to two differing noise depths: 1500 deg^ at the final 
SPT survey depth of 18 /iK-arcmin at 150 GHz, and 1000 deg^ at a depth of 54 /iK-arcmin. Clusters 
were selected on the basis of their SZ signal-to-noise ratio (S/N) in SPT maps, a quantity which has 
been demonstrated to correlate tightly with cluster mass. The S/N thresholds were chosen to achieve 
a comparable mass selection across survey fields of both depths. Cluster redshifts were obtained with 
optical and infrared imaging and spectroscopy from a variety of ground- and space-based facilities. 
The redshifts range from 0.098 < z < 1.132 with a median of Zmed = 0.40. The measured SZ S/N 
and redshifts lead to unbiased mass estimates ranging from 9.8 x 10^^ M^h^^ < M2oo(Pmean) < 
3.1 X 10^^ Mq hj^. Based on the SZ mass estimates, we find that none of the clusters are individually 
in significant tension with the ACDM cosmological model. We also test for evidence of non-Gaussianity 
based on the cluster sample and find the data show no preference for non- Gaussian perturbations. 
Subject headings: galaxies: clusters: individual, cosmology: observations 
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1. INTRODUCTION 

Galaxy clusters are the most massive collapsed objects 
in the Universe, with masses that range from 10^^ Mq 
to over 10^^ M©. Their abundance as a function of mass 
and reds hift can be used to constrain cosmological pa- 
rameters jWang fc Steinhardt| |1998^'Haiman et al.||2QQll 
Holder et al.||:^0011 mattyeTrW eller 2003; Mol nar et b1\ 
2004 r Wang et a H2004[|Lima fc H u 2007), and this con- 
straining power has now been demonstrat ed with real 
cluster samples i dentified in optical (e.g ., Rozo et al. 
[2010| ), X-ray (e.g., jVikhlinin et al.||2009b|), and, most re- 



cently, millim eter (mm) ( | Vanderlinde et al.||2010 Sehgal 



et 



Qtly, mnlir 
al.|2010 b") 



data. The most massive clusters are ot par 
ticular interest, especially at high redshifts. As tracers of 
the most extreme tails of the cosmological density field, 
these clusters can be used to place limits on the Gaus- 
siani ty of the initial density perturbations of the Universe 
(e.g., Matarrese et al.|2Q0Q[ ). Furthermore, massive, high- 
redshitt clusters provide laboratories for the study of as- 
trophysics (particularly galaxy formation and evolution) 
in dense environments in the early Universe. 

Although the largest existing catalogs of galaxy clus- 
ters are derived from optical and X-ray observations, 
clusters can also be identified by their interaction with 
cosmic microwave background (CMB) photons. The 
thermal Sunyaev-Zel'dovich (SZ) effect is a spectral dis- 
tortion of the CMB cause d by inverse- Compton s catter- 
ing with hot cluster gas ([Sunyaev fc Zel'dovich 1972). 
The surface brightness of tne ettect is independent oFred- 
shift, and the integrated thermal SZ effect from a cluster 
is expected to trace cluster mass with low scatte r (Bar- 

Carlstrom'"2001J [Motl eTaE] 



bosa et al."1996; Holder 



2005||Nagai e t al. 2007; Sta nek et a l. 2009), implying that 
cluster surveys should deliver nearly mass-limited cat- 
alogs of clusters to arbitrarily high redshift. With the re- 
cent development of bolometric receivers with hundreds 
or thousands of pixels, dedicated mm- wave SZ surveys 
over large areas of the sky are now bein g carried out by 
the South Pole Telescope (SPT, |Carlstrom et al ||2011|) 
and the Atacama Cosmology Telescope (ACT, lFowle£| 
eraI][2QQ7| ). Such surveys promise to be powerful tools 



tor cluster cosmology. 

The SPT is currently surveying 2500 deg^ of the south- 
ern sky at 95, 150, and 220 GHz. To date, roughly 
1500 deg^ have been observed to a depth of 18 /iK-arcmin 
at 150 GHz.^^ Motivated by preliminary evidence of 
surprisingly massive, high-redshift clusters in these first 
1500 deg^, we conducted "preview" observations of the 
remaining ~ 1000 deg^ of the SPT survey field during a 
three week period of the 2010 Austral winter, mapping 
this region to a noise level three times higher than the full 
survey depth (54 /iK-arcmin at 150 GHz). In this paper, 
we present a catalog of the 26 most significant galaxy 
clusters in the fuh 2500 deg^ SPT survey field and test 
whether the cluster masses and redshift distribution are 
consistent with those expected in a ACDM cosmology. 



timore, MD 21218 
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33 In this work, "/iK-arcmin" refers to the rms noise in equivalent 
CMB fluctuation temperature in a map with square 1' x 1' pixels. 



We complement our SZ cluster catalog with data from 
observations at other wavelengths. Spectroscopic or pho- 
tometric redshifts were obtained for each cluster as a part 
of a dedicated optical and infrared (IR) follow-up cam- 
paign. X-ray luminosities were also determined for each 
cluster using a combination of pointed observations with 
the Chandra satellite and measurements from the Roent- 



gensatellit (ROSAT) mission ( |Voges et al.||1999[ ) 



[2] de- 

peline, 



This paper is presented as follows: Section _ 
scribes the SPT observations, data reduction pipeline, 
and cluster-finding methodology. The catalog is pre- 
sented in Section [S] Section [4] introduces the simula- 
tions carried out to test purity and completeness, and to 
determine the scaling between the observable quantity 
(signal-to- noise ratio, S/N) and cluster mass. Section 
[S] describes the optical and IR follow-up measurements 
used to determine the cluster redshifts, and Section [6] 
presents X-ray luminosities for these clusters. Cosmo- 
logical implications of this cluster catalog are discussed 
in SectionlTl and our conclusions are presented in Section 

[HI 

Unless otherwise noted, we have assumed^ a 

WMAP7+BAO+i^o ACDM cosmology ( KomaTsu| 
"2QTT| with VLm = 0.272, = 0.728 and 

70.2 km s~^ Mpc~^ with distance measure- 
ments from Baryon Acoust ic Oscillations (BAO ) in 
the distribution of galaxies fPercival e t al.||2010) and 
" " (Ho) ■ ^ 



et al. 



^0 



the 



et al. 



Hu bble constant (i^o) measurement from |Riess 
( |2009[ ). Cluster mass estimates are reported 
irTTermsoi M2oo(Pmean), the mass enclosed within a 
radius corresponding to an average density of 200 times 
the mean density of the Universe. For the purposes 
of comparison with certain scaling relations in the 
literature, we also convert these masses into M5oo(Pcrit)5 
or the mass enclosed within a radius corresponding to 
an average density of 500 times the critical density. 
The conversion factor for each cluster is calculated 
assuming a Navarro- Frenk- Wh ite density profile a nd the 
mass-concentration relation of Duffy et al. (2008). 



2. INSTRUMENT, OBSERVATIONS, AND DATA 
REDUCTION 

2.1. The South Pole Telescope 

The SPT, a 10-meter off-axis Gregorian design with a 1 
deg^ field of view, has been searching for galaxy clusters 
in the mm- wave sky since its commissioning in 2007. The 
SPT is located within 1 km of the geographical South 
Pole. At an altitude of 2800 meters above sea level, the 
South Pole is one of the premier locations for mm-wave 
astronomy. The high altitude and low temperatures en- 
sure an atmosphere with low water-vapor content and 
excellent transparency. Meanwhile, the location near the 
Earth's rotational axis allows 24-hour access to the target 
fields. 

The SZ receiver currently mounted on the teles cope 
consists of 960 transition-edge-sensor bolometers (Lee 
et al. 1998[), cooled to a temperature of 280 mK. These 
bolorneters are split into six wedges each containing 160 
detectors. The sensitivity and configuration of these 
wedges have changed over the four years of scientific op- 
eration. In 2007, we fielded a preliminary array with 
wedges at all three frequencies but with limited sensitiv- 
ity. In 2008, the array contained a single 95 GHz wedge, 
three 150 GHz wedges and two 220 GHz wedges. The 
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95 GHz wedge did not produce science-quality data, but 
the 150 and 220 GHz wedges performed to specification. 
In 2009, the 95 GHz wedge was replaced with a wedge 
with much higher sensitivity, and one of the 220 GHz 
wedges was replaced by a 150 GHz wedge, resulting in 
an array with one wedge at 95 GHz, four wedges at 150 
GHz, and one wedge at 220 GHz. The focal plane config- 
uration has remained the same since 2009. The 10-meter 
primary is conservatively illuminated, resulting in beam 
sizes (FWHM) of approximately 1.6', 1.1', and 1.0' at 95, 
150, and 220 GHz. 

The SPT team has previously published tw-o SZ- 
selected cluster samples: Staniszewski et al. (2009 here- 
after S09) presented four clusters (including three newly 
discovered clusters) selected from ^ 40 deg^ of 200 7 and 
2008 150 GHz data, while Vanderlinde et aH 020101 here- 
after VIO) presented 21 clusters (including 12 new dis- 
coveries beyond S09) selected from 150 GHz data in the 
full 200 deg2 of 2008 observations. 

2.2. Observations 

The SPT-SZ survey area is a contiguous 2500 deg^ re- 
gion defined by the boundaries 20h < R.A. < 24h; Oh < 
R.A. < 7h and —65° < S < —40°. This region comprises 
the majority of the low-dust-emission southern sky be- 
low 6 = —40°. Observing fields north of this declination 
becomes difficult with the SPT because of the increased 
atmospheric loading and attenuation at low elevation. 
We split the survey region into 19 fields, ranging in size 
from ~ 70 deg^ to ~ 230 deg^. Single observations of 
fields of this size can be completed in an hour or two, 
allowing a regular sc hedule of inter leaved calibrations 
(see S09 and Carlstrom et al.| ( |201lD for details). The 
exact size, shape, location, and order of observation of 
these fields are determined by a combination of factors 
including availability of data at other wavelengths, sun 
avoidance (some of our observations take place during 
the Austral spring and summer), and the desire to have 
a final survey area that is easy to define. The location, 
year of observation and size for each field are shown in 
Table [l] Figure [l] shows the full survey region and the 
individua l field borders ov erlaid on the 100 jam dust map 
from Schlegel et aT] ( [199 81). 

The standard operating mode of the SPT is to ob- 
serve a target field by scanning back and forth in az- 
imuth across the field followed by a step in elevation. 
These steps are large compared to the beam size, so sub- 
sequent observations of the field have a small offset in 
elevation applied in order to oversample the sky. Cer- 
tain fields were observed in what is called a "lead-trail" 
mode. In this observing mode, the lead half of a field 
was scanned followed by the trail half, as opposed to 
scanning the entire azimuth range of the field in a single 
scan. This strategy was employed to safeguard against 
possible ground contamination, but we see no evidence 
of such contamination on the angular scales of interest to 
this work. Approximately two-thirds of the ra21hdec-50 
observations^^ were taken using an elevation scan mode 
rather than scanning in azimuth. In this mode, the tele- 
scope was parked at a fixed azimuth and scanned up and 
down in elevation, allowing the sky to drift through the 

Coordinates in the field names refer to the R.A. and d of the 
center of the field. 
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Fig. 1. — Outlines of the SPT-SZ survey fields overlaid on an or- 
thographic projection of the IRAS 100 fim dust map from |Schlegel| 
et al. (1998). The sky is rotated such that the South Celestial Pole 
is at the top of the globe, and R.A. = Ih faces the viewer. Green 
lines indicate fields observed in 2008, red lines indicate fields ob- 
served in 2009, blue lines indicate fields observed to full depth in 
2010, and yellow lines indicate fields observed to preview depth in 
2010, which will be completed to full depth in 2011. 

field of view. We include data from both azimuth and 
elevation scans on this field. We have investigated the 
effects of these different scan strategies on noise prop- 
erties and cluster finding and found these effects to be 
negligible. 

As mentioned in the introduction, observations have 
been completed to full survey depth (a noise level of 
18 /iK-arcmin at 150 GHz) for roughly 1500 deg^ of the 
SPT survey region. Initial preview observations of the 
remaining ~ 1000 deg^ were performed in late 2010, to 
a noise level of 54 /iK-arcmin at 150 GHz, or three times 
the full-depth noise level. The results presented in this 
work are based on data from both the full-depth and 
preview-depth fields. 

One field, the ra23h30dec-55 field, was observed in 
both 2008 and 2010. Given the higher quality of the 
95 GHz wedge in 2010, that year's data is us ed in pref- 
erence to the 2008 data, which was used in V10[ For 
this reason, the properties of S PT-C L J2337-5942 are 
not identical to those reported in VIO We choose not to 
combine the 2008 and 2010 data in order to reduce the 
number of different map depths considered. 



2.3. Data processing, calibration, and map-making 

The data reduction pipeline applied to the SPT data 
is very similar to that des cribed in previous SPT pa- 
pers such as S09, VIO, and [Shirokoff et al.| (12011"). An 
overview of the processing is presented here, highlighting 
differences with earlier SPT releases. The same data- 
processing and map-making procedure is used for each 
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TABLE 1 
2008-2011 SPT Fields 



Field Name 


Obs Year 


Area [deg^] 


raondUdec-oo 


2008 


90 


radnc5Udec-oU 


2009 


230 


razlndec-5U 


2009 


200 


razlndec-oU 


2009 


150 


raUnoUdec-oU 


2010 


160 


ralndec-oLl 


2010 


150 


raznoUdec-5U 


2010 


160 


ra4hl0dec-50 


2010 


160 


ra5h30dec-45 


2010 


110 


ra23h30dec-55 


2010* 


100 


ralhdec-42.5 


2010S 


110 


ra3h30dec-42.5 


2010S 


170 


ra6hdec-62.5 


2010S 


70 


ra6h30dec-45 


2010S 


110 


ra6h30dec-55 


2010S 


90 


ra21hdec-42.5 


2010S 


110 


ra22h30dec-55 


2010S 


80 


ra23hdec-45 


2010S 


210 


ra23hdec-62.5 


2010S 


70 



Note. — The field centers, year of ob- 
servation and area of the SPT fields. The 
19 SPT fields cover a total of roughly 2500 
deg^. The nominal noise level of an SPT 
field is 18 /iK-arcmin. "2010S" refers to 
those fields observed to three times the 
survey noise level (54 /iK-arcmin). The 
ra23h30dec-55 field was also observed in 
2008, but only the 2010 data is used in this 
work. 

field, with minor adaptations in filtering to produce uni- 
form map properties regardless of scan strategy. 

The first steps in processing are to flag regions of com- 
promised data (for instance, time samples with cosmic 
ray events) and to reconstruct the pointing for each de- 
tector. We then calibrate the time-ordered data (TOD) 
to CMB temperature units. As in S09, this calibration is 
based on observations of a galactic HII region (RCW38) . 
The TOD is filtered and co-added into the final single- 
frequency map with inverse-noise weighting. 

The filtering consists of bandpass filtering the TOD 
and removing correlated noise between detectors. The 
high-pass filter is implemented by removing a ninth-order 
Legendre polynomial and a set of Fourier modes from 
each scan. The highest-frequency Fourier modes removed 
correspond to an angular frequency of k = 400 in the 
scan direction. Depending on the scan strategy used 
for the observation, this filter acts as a high-pass filt er in 
the R.A. or decl. direction. This differs slightly from |VlQ 
where only a first-order polynomial was removed, and the 
set of Fourier modes removed was defined by temporal 
frequency (/ < .25 Hz) rather than angular frequency 
{k < 400). The cutoff definition was altered to handle 
variable scan speeds. For the 2008 scan speeds, 0.25 Hz 
corresponds to k 360, which means that the /c-space 
high-pass cutoff is slightly higher in this work than in 
VIO A low-pass filter was also applied (with a cutoff 
at k 30000) to avoid aliasing of high-frequency TOD 
noise when the data is binned into a map. 

Atmospheric noise is correlated across the entire focal 
plane. VIO removed the mean and slope across all detec- 
tors in a frequency band at each time sample. However, 

We use the fiat-sky approximation throughout this work, so 
|k|^l 



the number of detectors at the two frequencies used in 
this work (95 and 150 GHz) differ by a factor of four, so 
this scheme would filter different spatial mod es on the 
sky a t each frequency. Instead, as was done in [Shirok offl 
'etarf20ri), the mean of the TOD across a geometrically 
compact set of one quarter of the 150 GHz detectors or all 
the 95 GHz detectors (i.e., across one detector wedge) is 
subtracted at each time sample. This acts as an isotropic 
high-pass filter with a cutoff at roughly k = 500. 

2.4. Cluster Finding 

As discussed in Section[2]2j most of the SPT fields have 
been observed in three frequency bands, centered at 95, 
150, and 220 GHz. Multiple sky signals and sources of 
noise contribute to each single-frequency co-added map 
of a field, and each of these contributions has unique spa- 
tial and spectral properties. Primary CMB fluctuations, 
emissive point sources, and noise (both atmospheric and 
instrumental) contribute to the maps at all three frequen- 
cies. A small signal from the kinetic SZ (kSZ) effect (due 
to the interaction between CMB photons and free elec- 
trons with a bulk velocity) also contributes to all three 
frequencies. Most importantly for this work, the 95 GHz 
and 150 GHz maps contain an additional signal due to 
the thermal SZ (tSZ) effect from clusters. Because we 
can predict the spectral signature of the tSZ effect (up 
to a small relativistic correction), we can combine the 
maps from the three bands to maximize sensitivity to 
tSZ and minimize noise and other contaminants. Fur- 
thermore, we can use the fact that the galaxy clusters 
we expect to find have a different spatial profile than 
other signals and noise to construct a spatial filter that 
maximizes sensi tivity to clus ter-shaped signals. 

As shown by Melin et al.| ([2006) and others, the op- 
timal^^ way to extract a cluster-shaped tSZ signal from 
our data is to construct a simultaneous spatial-spectral 
filter. We begin by assuming that the maps are fully 
described by 

T(x,^,)= (1) 

5(X, Ui) * [fsz{^i)TcMBysz{^) + nastro(x, l^i)] 
H~ '^noise(x, i^i), 

where ysz is the true tSZ sky signal in units of the Comp- 
ton y parameter, Tcmb is the mean temperature of the 
CMB, /sz encodes the frequency scaling of the tSZ effect 
relative to primary CMB fluctuations (e.g., Ca rlstrom| 
|et al.||2002[ ), nastro and nnoise are the astrophysical sig- 
nals and instrument/atmospheric noise we wish to de- 
weight, B(^^Vi) encodes the instrument beam and any 
filtering applied in the analysis, and denotes con- 
volution. Given this assumption, the matched spatial- 
spectral filter is given by 

%l^{k:j,,ky,Vi) = cr~'^ ^'^^^{kx^ky)fsz{i^j)Smt{kx,ky,Uj). 
j 

(2) 

This method is in fact optimal only under certain assump- 
tions, the most important of which are that all sources of noise 
and unwanted astrophysical signals are random and translation- 
ally invariant, and that the exact spectral and spatial behavior of 
every component of signal and noise are known perfectly. 
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Here, is the predicted variance in the filtered map 
cr^^ = ^fsz{j^i)Smt{kx^ ky, Vi) ^~^(k:,, ky) X (3) 

fsz{^j)Smt{kx,ky,Uj), 

and ^fiit is the assumed cluster profile convolved with 
B((K^Vi). The kx and ky arguments are included explic- 
itly in ^fiit because, while the underlying cluster profile 
is assumed to be a zimu thally symmetric, the filtering de- 
scribed in Section [23| is anisotropic. The v argument is 
included explicitly to account for the fact that the filter- 
ing and instrument beam can be different in the different 
SPT bands. N is the band-band, pixel-pixel covariance 
matrix describing the noise and non-tSZ signal 

Nabij = (4) 
^astro ^noise 
[5(X5, Vj) * nastro(x5, ^j) + nnoise(x6, ^j)])- 

Under the assumption that these components are 
translationally invariant, the pixel-pixel part of this ma- 
trix will be diagonal in the Fourier domain, which is why 
we only include the band indices in Equation [2] This 
also means that '0 can be evaluated separately at each 
value of {kx^ky}^ and the largest matrix that needs to 
be inverted is A^bands-by-ATbands. 

There should be no correlation between the astrophysi- 
cal signals and instrumental/atmospheric noise, in which 
case N can be separated into Nastro and Nnoise- Further- 
more, the instrumental noise should be uncorrelated be- 
tween bands, although the atmospheric noise may have 
correlations. We have performed correlation analyses on 
SPT maps similar to those used in this work and found 
little, if any, noise correlation between bands, which is 
expected because the correlated part of the atmospheric 
emissio n is largely removed in the filtering described in 
Section 2.3 In this case, we can estimate Nnoise individu- 
ally in each b and. We do so using the jackknife procedure 
described in,VTO|and S09. 

Our model for Nastro is a combination of primary and 
lensed CMB fluctuations, point sources below the SPT 
detection threshold, kSZ, and tSZ from clusters below the 
SPT detection threshold. The power-spectrum shapes 
and 150 GHz amplit udes f or these components are iden- 
tical to those used in |V1Q[ The spectral behavior of the 
primary CMB, kSZ, and tSZ components are known (up 
to the relativistic correction for tSZ, which we ignore). 
The spectral behavior of the point sources is assumed to 
be such that the flux density of a given source follows a 
power law in frequency {v/vq)'^^ with a = 3.6. This is 
consistent with the behavior of dusty, star-forming galax- 
ies (DSFGs) below the SPT detection threshold. Radio 
sources below the SPT detection threshold are expected 
to contri bute ne gligibl y to the map rms compared to 



noise ( |Hall et al.|m jShirokoff et al.||2Qll| . 

As in |VlU| and 5UMhe source template S is described 
by a projected spherical /3-model, with f3 fixed to 1, 



AT = ATo(l + 6>V6>c)" 



(5) 



where the normalization ATq and the core radius Oc are 
free parameters. As in VIO and S09, twelve different 
matched filters were constructed and applied to the data. 



each with a different core radius, spaced evenly between 
0.25' and 3.0^ As in V1Q[ point sources detected above 
5(7 were masked out to a radius of 4', with the value 
inside that radius set to the average of the surrounding 
pixels. For extended sources, a custom mask was applied, 
covering the shape of the emission. In both cases, cluster 
detections inside the masked region or within 4' of the 
outer edge of the masked region were rejected, meaning 
that any detection within 8' of a masked point source 
was rejected. 

The method of extracting clusters from the filtered 
maps in this work (including S/N estimation and peak 
detection) is identical to that used in ,V1Q,, As in VIO, 
we refer to the detection significance maximized across 
all twelve matched filters as ^, and we use ^ as the pri- 
mary SZ observable. 

We have measured the performance of the multi-band 
cluster detection algorithm relative to single-band detec- 
tion and found significant improvement when we add the 
95 GHz data — which has roughly a factor of two higher 
noise in CMB fluctuation temperature than the 150 GHz 
data — over 150 GHz data alone. We find very little fur- 
ther improvement when we add the 220 GHz data, which 
has a factor of five higher noise than the 150 GHz data. 
As the 220 GHz SPT maps are not currently deep enough 
to add significantly to the cluster detection efficiency, we 
use only 95 and 150 GHz data for all the results presented 
in this work. 

For one field, ra5h30dec-55, only 150 and 220 GHz data 
exist, and we only perform single-frequency 150 GHz 
cluster finding on this field's map (with results indis- 
tinguishable from those reported in |V10| ). Using the 
measured improvement in cluster S/N in multi-band vs. 
single-band data (roughly 20% averaged across all red- 
shifts and cluster sizes), we determined that the most 
significant cluster in this field would not have made it 
into the catalog in this work even if we had 95 GHz data 
on this field (see Section [s] for details of the selection of 
clusters for this work). 

3. CATALOG 

The cluster-finding pipeline described above returns 
hundreds or even thousands of candidates (depending on 
the threshold value of ^) within the 2500 deg^ survey 
area. Some of these candidates have already been re- 
ported in S09 and VIO, and upcoming publications will 
continue to expand the SPT catalog, including candi- 
dates down to S/N values as low as ^ = 4.5, where the 
purity of the sample is still estimated to be well over 
50%. The aim for this paper, however, is to search the 
full survey area for the clusters that have the greatest po- 
tential to test the current cosmological model — in other 
words, the most massive clusters. At any given redshift, 
the most massive clusters will correspond to the clusters 
with the highest SZ significance, so we present here a 
catalog of the 26 most significant detections in the 2500 
deg^ survey area. 

This catalog is constructed by setting a high signif- 
icance threshold of ^ > 7 for all shallow fields. This 
threshold was chosen to ensure zero false detections at 
high confidence (see Section |4] for details on the false 
detection rate estimate) and to limit the scattering of 
low-mass systems into the sample. To choose a deep- 
field threshold to match this shallow-field threshold, we 
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use simulated observations (described in Section H| at 
both field depths using the same underlying simulated 
SZ skies. For each cluster which met the ^ > 7 thresh- 
old in the shallow-field simulations, a deep-field ^ was 
calculated. The shallow-to-deep ^ ratio was then calcu- 
lated for each cluster and found to vary from 0.4-0.6, 
depending on cluster characteristics. These ratios were 
averaged to yield an approximately equivalent threshold 
value of ^ > 13 for the full-depth fields. Note, how- 
ever, that significance values do not scale simply with 
field depth. The matched-filter noise is a combination 
of the CMB, point sources and observational noise, and 
the relative contribution of each is /c-dependent, varying 
significantly between the different depths. For example, 
in shallow fields the relative contribution of the CMB is 
less, tilting the matched filter to prefer larger scales; cor- 
respondingly, it prefers smaller scales in deep fields. A 
simple, direct translation of "preview-depth ^" to "full- 
depth ^" is therefore not possible on a cluster-by-cluster 
basis; the matching of thresholds is approximate and de- 
signed only to match the average scaling. 

Nevertheless, as shown in Figure |2j the rough 50% 
and 90% completeness contours using the two depths are 
fairly well matched, particularly at high redshift. These 
contours are calculated as in VIO, by matching detec- 
tions in simulated observations with halos identified in 
the underlying dark matter simulations in bins of mass 
and redshift. For the high ^ thresholds in this work — 
which correspond to very massive dark matter halos — 
the number of matching halos is small in most bins, so 
the completeness contours are poorly sampled and noisy, 
but the general agreement between the two sets of con- 
tours is clear. Note that the uncertainty in these curves 
does not directly affect the remainder of this work, as the 
selection function is always considered as a strict thresh- 
old in ^; the selection as a function of mass and redshift 
is presented in Figure [2] purely for illustrative purposes. 

For each cluster, both the redshift, and the X-ray lumi- 
nosity were determined (as outlined in Section [sl and Sec- 
tion |6|. The catalog is presented in Table [6] and thumb- 
nail images of each cluster in the SZ and optical/IR are 
shown in Appendix A. 

3.1. Notable Clusters 

Of the 26 clusters reported in this paper, 12 are new 
discoveries, one was previously reported in VIO, and 13 
others have been identified in other optical. X-ray, and 
SZ cluster catalogs, with 7 having multiple identifica- 
tions. The previously identified clust ers include seven 



clusters in the optical Abell catalog (Abell et al. 1989 



nine clusters in the X-ray ROSAT-ESO t'lux Limited 
ray Galaxy c luster survey catalog (REFLEX, Bohringer 



et al. 2004a), and six clusters in the mm-wave Ata- 
cama Cosmology Telescope catalog (ACT, [Marriage etal. 



|20TT1). These cross-associations and alternative identiti 
cations are noted in Table [6j In this section, we discuss 
particularly notable clusters in the SPT catalog. 
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Fig. 2. — Completeness contours as a function of mass 
(M2oo(p mean)) and rcdshift for both sets of map depths and ^ 
thresholds, estimated from simulated observations. 50% and 90% 
contours are shown for shallow fields with ^ > 7 (black, solid lines) 
and for deep fields with ^ > 13 (red, dashed lines). 

is completed to full depth. Given the redshift of this 
cluster {z = 0.78), it is expected to be one of the rarest 
objects in the SPT survey (see Figure [5|. 

SPT-CL J 06 15-5746— This cluster has the second high- 
est redshift of any cluster in this paper, with a redshift 
oi z = 0.972. Based on its ROSAT faint source catalog 
counterpart, it is measured to be the fourth most X-ray 
luminous cluster in this catalog (see Section [6|. 

SPT-CL J0658-5556— This cluster is the well-known 
Bullet cluster, otherwise known as lES 0657-558. It has 
been extensively s tudied in multiple wavelengths (e.g., 
Clowe et al. 2006) and is known to be one of the most 
massive and X-ray-luminous clusters in the Universe. It 
is expected to be the most massive cluster in the final 
SPT catalog. 

SPT-CL J2 106-5844— Multi- wavelength observations of 
this SPT-discove red cluster are discussed in detail in lFo-l 
ley et al. ( |2011 ). This is the highest-redshift cluster 
[z = 1.132) spectroscopically confirmed in the SPT sur- 
vey. X-ray observations from Chandra measure an X -ray 
luminosity of [0.5-2.0 keV]= 13.9 x 10"^"^ erg s-^ 



ley et al. 2011[ ), comparable to the X-ray luminosity of 
the Bullet cluster. The mass we report in Table |6] for 
SPT-CL J2 106-5844 is slightly (5-10%) d iscrepant with 



the SZ mass reported in Foley et al.| (2011 ), although the 
difference is much smaller than either value's la uncer- 



tainty. The difference arises because ""Foley et al.J ([2011 ) 
use the single-band 150 GHz ^ and the exact VIO scaling 
relation to derive the mass, whereas this work uses the 
multi-band ^ and the multi-band scaling re latio n devel- 
oped specifically for this work. (See Section 



4.1 



and VIO 



SPT-CL J0102-m^ This cluster was first reported in ^^^^-^^ ^^^^ estimation and scalingVelations 



Marriage et al. (2011 ). It is the most significant detection 
to date in the full SPT survey by nearly a factor of two. 
It has a comparable X-ray luminosity and beam-averaged 
SZ decrement to the Bullet cluster and AS 1063, whose 
SZ significances should be similar when the SPT survey 



SPT-CL J2248-4431— This cluster is also known as 
AS1063. It is the second most X-ray luminous cluste r 
in the REFLEX X-ray survey ("Bohringer e t al.||2004a] , 
even more luminous than the Bullet cluster. It has the 
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second highest estimated mass for any cluster in this pa- 
per. 

SPT-CL J2344-4243— From its redshift and its ROSAT 
bright source catalog counterpart (see Section [6|, this 
cluster is measured to have the largest X-ray lummosity 
of any cluster in this paper. A bright Type 2 Seyfert 
galaxy at redshift 0.5975, 2MASX J23444387-4243124, 
is located 19'^ from the SZ cluster centroid. This redshift 
is consistent with our photometric red-sequence redshift 
estimate of 0.62 for the cluster (see SectionT5|, suggesting 
this galaxy may be in or near this cluster. 

3.2. Point Source Veto 

As discussed in Section [2^ emissive sources above 5cr 
are masked in the cluster finding procedure, and any clus- 
ter detections within 8' of a masked point source are re- 
jected, because residual source fiux or artifacts due to the 
masking can cause spurious decrements when the maps 
are filtered. However, this rather conservative procedure 
can result in rejecting a cluster detection that was only 
marginally affected by the nearby emissive source. 

To test this scenario, we re-ran the cluster finding algo- 
rithm on all the fields used in this work with only the very 
brightest (6'i5oghz > 50 mJy) sources masked (as com- 
pared to the 5(7 thresholds of roughly 6.4 and 17 mJy in 
the full-depth and shallow fields). Each detection above 
the ^ threshold for this paper was visually inspected, and 
the vast majority were rejected as obvious point-source- 
related artifacts. However, two objects were clearly real 
detections. One of these detections, with ^ = 19.3 (in a 
full-depth field), is a known cluster (Abell S0295), and we 
include it in our catalog as SPT-CL J0245-5302. There 
are two > 5cr sources within 8 arcmin of this cluster, 
but neither is strong enough to affect the ^ measurement 
by more than 1 or 2(7. Furthermore, the sources lie to 
the north and south of the cluster, and wings from the 
matched filter are predominantly along the scan direction 
(east- west in all SPT fields but one). However, because 
this cluster was not found by the original version of the 
cluster finding algorithm, we do not include it any cos- 
mological analysis. 

The other detection, ^ = 13.3 (full-depth), has an Sa 
source almost directly to the east. Because this cluster 
(SPT-CL J2142-6419, which wih likely appear in a future 
catalog) could have been bumped over our ^ = 13 full- 
depth threshold by the filtering wing of the point source, 
we choose not to include it in this work. 

4. SIMULATIONS 

Simulated observations were used to characterize the 
catalog presented in this work. These simulations were 
modeled after those used in |V10| where full details can 
be found; a brief summary is provided here. 

We create random Gaussian realizations of the power 
spectra of primary CMB anisotropy, kSZ, and point 
sources below the SPT detection threshold. The pri- 
mary CMB power spectrum is chosen to match the 
WMAP7 be st-fit ACDM model (L arson et al. 2011) , and 
we add the Sehgal et al.| ( 2010a| predicted "kSZ power 
spectrum to the expected CMB anisotropy. The point- 
source power spectrum includes terms corresponding to 
synchrotron-dominated sour ces and DSFGs with a mpli- 
tudes based on the results of Shirokoff et al. (2011). We 



include Gaussian realizations of the power from Poisson 
distributions of both sources with amplitudes at 150 GHz 
and k = 3000 of 1^3000 — 1-3 /iK^ for the synchrotron 
sources and D^^^^ = 7.7 /iK^ for the DSFGs. We as- 
sume a spectral index of a = —0.6 for the synchrotron 
sources and a = 3.6 for DSFGs. In addition to the Pois- 
son power, we model a clustered DSFG component with 
an angular multipole dependence of oc k and ampli- 
tude at 150 GHz and k = 3000 of D^qoo = 5.9 /iK^. SZ 
sk ies were sim ulated as in VIO] using the methodology 
of Shaw et al.| (^2009); the fiducial simulations from that 
work are used here. Briefiy, they consist of semi-analytic 
gas models pasted over halos identified in N-body dark 
matter simulations. At each frequency (95 and 150 GHz), 
forty 100 deg^ sky maps were generated. 

These simulated skies were processed with an analytic 
approximation to the SPT transfer function, consisting 
of a Gaussian beam, an isotropic high-pass filter and a 
high-pass filter along the R.A. direction. The filters were 
arranged to model t he eff ect of the SPT data processing 
described in Section [231 with the R.A. high-pass set to 
match the /cr.a. = 400 cutoff of the Fourier mode removal 
and the isotropic high-pass set at /c = 500 to approximate 
the spatial template removal. 

Noise realizations were generated at each frequency for 
both the full survey and preview depths. The noise power 
was measured in differenced (jackknife) maps for each 
field, and these powers were averaged across the set of 
fields at each frequency and depth. Forty Gaussian ran- 
dom realizations of each of these averages were then gen- 
erated and added to the processed simulated sky maps. 

These simulations were then subjected to the same 
cluster-finding pipeline applied to the real data, and re- 
covered clusters were matched with the underlying cata- 
log of massive halos associated with the SZ simulations. 

4.1. Mass Scaling Relations and Unbiased Mass 
Estimates 

A number of different techniques are available for ob- 
taining cluster mass estimates from SZ measurements. 
The integrated SZ fi ux, Y, is expected to be a ti^ht proxy 
for the cluster mass (^Barbosa et al. 1996 Holder & Carl- 
strQm|[ 2QQr; 'M otl et al IJl^OOSI [Nagal et al.||l^007| l^haw 
e^al.| |2008; Stanek et al!][2QQ9|. Unfortunately, the dlSi^ 
culty in determining the correct filter scale Oc from SPT 
data alone adds significant scatter to the scaling relation 
of mass with Y (see VIO for details). In simulations, 
the SPT significance ^ has a smaller scatter than our cur- 
rent integrated Y estimates, and, as in VIO, we use the 
significance as a mass proxy. 

Due to the significant impact of noise biases, a direct 
^-M scaling relation is complex and difficult t o cha rac- 
terize. Instead, following the prescription of |V10[ we 
introduce an unbiased significance whose scaling with 
mass M2oo(Pmean) takcs the form 



M 



5 X W^MqH-^ 



1.6 



(6) 



where A is a normalization, B a mass evolution and C a 

Effort is currently underway to use our mult i- frequency data 
to improve the determination of 6c- 
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TABLE 2 

Mass Scaling Relations 



Depth 


A 


B 


C 


Scatter 


Full survey 


7.50 


1.32 


1.64 


0.21 


Preview 


3.50 


1.29 


0.87 


0.16 



redshift evolution. In simulated maps of both depths, 
( was calculated for each cluster as in VI 0, by deter- 
mining the preferred filter scale and cluster position in 
the absence of noise then averaging the detection signif- 
icance at that filter scale and position over many noise 
realizations. Mass scaling relations were fit to the subset 
of these with M > 2 x lO^'^/i-^M© and z > 0.3 by min- 
imizing the residual logarithmic scatter in C about the 
relation. These relations are given in Table|2l As these 
are based on the same SZ simulations used in VIO, they 
can be viewed as equivalent to the relations presented in 
that work. 

Uncertainties in the SZ modeling lead to significant 
systematic uncertainties on these scaling relation param- 
eters. Following V1Q[ we apply conservative 30%, 20%, 
50%, 20% Gaussian uncertainties to A, B, C, and scatter, 
respectively. 

Mass estimates are constructed as in |V10| with slight 
modifications to account for the different field depths. 
Details of the conversion from ^ to mass are given in 
Appendices B and D of VIO. Briefly, we calculate the 
conditional probability of detecting a cluster of mass M 
at a given value of ^, P(^|M), and then apply a mass- 
function prior to create the posterior probability P(M|^). 
This procedure accounts for two types of bias in the mass 
estimate, the first due to the fact that we have maximized 
^ over many filter choices and positions in the map, and 
the second due to the combination of the steepness of the 
cluster mass function and observational noise or scatter 
in the mass-observable scaling-relation. The latter ef- 
fect, which results in more low-mass systems scattering 
up into a given ^ bin than high-mass systems scatter- 
ing down, is related to the phenomenon of Eddington 
bias'^o (Eddington 1913). For the very high detection 
significances used in this work, the maximization bias is 
completely negligible compared to the bia s due to this 
asymm etric scatter. As in VIO, we use the [Tinker et al. 
([2008^ mass function evaluated at the maximum likeli- 
hood point in the WMAP7+V10 chain as our prior. This 
method produces unbiased posterior estimates for cluster 
masses, assuming the validity of the simulations and of 
the various priors app lied. 



In Andersson et aL\ (2010), we compared SZ inferred 
masses calculated witn this method to X-ray mass esti- 
mates for the 15 clusters from VIO that had X-ray mea- 
surements. Overall, we found agreement between the SZ 
and X-ray mass estimates near the quoted level of the 

Note that, for consistency with VIO, this relation is given in 
terms of h (i.e., /iioo), not hro 

The redshift cutoff is due to the fact that it was found in 
VIO that the power-law parametrization of the scaling relation fails 
to fully capture the behavior of the SPT selection function below 
z = 0.3. 

Strictly speaking, "Eddington bias" refers to the bias in num- 
ber counts caused by this asymmetric scatter. The bias in the 
measured properties of individual objects is sometimes erroneously 
referred to as "Malmquist bias" . 



systematic uncertainties of the SZ mass estimates. How- 
ever, there was a significant statistical offset, with the 
SZ-inferred masses lower by a factor of 0.78 ± 0.06 aver- 
aged over the sample, which we do not correct for in the 
masses in Table IH This factor could have a redshift or 
mass dependence and naively applying a correction fac- 
tor would ignore these effects. There was some evidence 
for this in , Andersson et al.J 2010), where the lowest red- 
shift and most massive cluster had the most discrepant 
SZ and X-ray inferred masses. We are currently pursu- 
ing an analysis that jointly constrains the SZ and X-ray 
mass- observable relations with cosmology ( [Benson et al^ 
2011), which will more accurately quantify any system- 
atic offset between the SZ and X-ray mass estimates. 
For this work, we consider the quoted uncertainty of the 
SZ mass estimates to be a reasonable estimate of the sys- 
tematic uncertainty, and note that there is some evidence 
that the SZ mass estimates are low by ~25%. We show 
in Section that none of the conclusions in our cosmo- 
logical analyses would change if we naively applied this 
scaling factor. 

Finally, we note that, although the scaling relation 
fits were only performed on simulated clusters above 
z = 0.3, we nevertheless report SZ-derived masses for 
several z < 0.3 clusters in Table [6l These mass esti- 
mates are extrapolations of the scaling relations to areas 
of parameter space in which they have not been tied to 
simulations and may therefore be subject to further sys- 
tematic uncertainties. For this reason, we do not use any 
clusters at z < 0.3 in the likelihood calculations described 
in Section [7^21 



4.2. Purity and Completeness 

To test the likelihood of false detections, simulated ob- 
servations were generated omitting the SZ signal, and run 
through the cluster-finding pipeline. The false detection 
rate was found to be a rapidly falling function of the de- 
tection threshold. No false detections were found above a 
significance of ^ = 6 in simulations of 4000 deg^ at either 
depth. Given that the lowest threshold used in generat- 
ing the catalog presented in this work is ^ > 7, it is highly 
improbable that it contains any false detections. This 
is confirmed by the multi-band followup, which shows 
counterparts for each cluster in the catalog. 

The catalog is complete above threshold ^ values by 
construction. As discussed in Section [3j two factors make 
it difficult to quantitatively convert this into a mass and 
redshift completeness. The clusters in this catalog lie at 
the extreme high-mass end, and, as such, are rare in the 
simulated skies, yielding insufficient statistics to obtain 
a robust estimate of their detect ability. Furthermore, 
there is large uncertainty associated with modeling the 
gas attached to halos in the simulation that makes any 
threshold uncertain at the ~ 30% level. Figure [2] shows 
our best estimate of 50% and 90% completenessTor the 
two sets of field depths and ^ thresholds. We note that, 
as the same simulated SZ skies were used for both depths, 
the variance due to the limited sample size will appear as 
a coherent shift, rather than a scatter, between the two 
sets of curves 

5. OPTICAL AND INFRARED DATA 

Multi-band imaging from both ground- and space- 
based facilities has been obtained for clusters in the cat- 
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Fig. 3. — Redshift histogram of the sample. 



alog, for the purpose of cluster redshift estimation where 
no previous spectroscopic redshift could be determined 
from the literature. We have also carried out new multi- 
slit spectroscopy on some of the clusters. A summary of 
the cluster redshifts is shown in Table [Sj and a descrip- 
tion of our methods is outlined below. Figure [3] shows 
the redshift histogram of our cluster sample. 

5.1. Ground-based Imaging 

The previously unknown clusters in this catalog were 
imaged with the cameras shown in Table The bwope 
1 m at Las Campanas Observatory, equipped with the 
SITe3 CCD detector and BVRI filters, provided sensi- 
tivity to clusters at z < 0.7. We chose at least two pass- 
bands that spanned the 4000 A break, as determined ini- 
tially using cluster redshifts estimated by eye from Digi- 
tized Sky Survey (DSS) images. We required a detection 
of 0.4L* early- type cluster galaxies at about 5<j. We 
then iterated on this strategy: if the cluster was still not 
sufficiently detected in the initial set of exposures, we 
updated our best redshift estimate using the new data, 
if possible, and reobserved. 

The Blanco 4 m at Cerro Tololo Inter- American Obser- 
vatory and the Magellan 6.5 m telescopes at Las Cam- 
panas Observatory, provided sensitivity to clusters at 
z <l. The Blanco/MOSAIC-II, Magellan/LDSS3, and 
Magellan/IMACS CCD cameras were used with griz fil- 
ters. The same iterative strategy was implemented until 
we reached our r equired d etection. 

As detailed in Foley et al.^ (2011J , 50 min of preimag- 
ing with the VLT's~FUR5z camera in / provided addi- 
tional broadband data for our highest-redshift cluster, 
SPT-CL J2106-5844, in our corresponding spectroscopic 
program. 

All images were reduced in a uniform manner using 
the same software and methods described in detail in the 
previous, closely related work of 



High et al.|(|2010] 



here- 
tellar 



the same red-sequence software described inlHlOl These 
redshifts are estimated to be accurate to (Jz/ z) ^ 2- 
3% (statistical plus sytematic), as determined using a 
larger subset of clusters with added spectroscopic red- 
shift data. An exception to this are the SWOPE derived 
redshifts (SPT-CL J0245-5302 and SPT-CL J041 1-4819) 
which are accurate to cFz/(X + ^) ^ 4-5%. This was due 
to using Johnson filters (as opposed to Gunn-Sloan fil- 
ters), and calibrating the photometry with a synthetic 
stella r locus (rathe r than the median SDSS stellar locus 
from [Covey et al"! ([2007|) using the PHOENIX stehar 
model atmosphere library (Brott & Hauschildt 2005). 

A parallel reduction of the Blanco/MOSAIC-II data 
was performed using th e Dark Energy Survey data man - 
agement system (jNgeow et al."2006 Mohr et al. 2008). 
Redshifts for theT2 clusters with MO!SAIC-II imagmg 
were independently measured from these data using Ar- 
tificial Neural Network (ANN^ Collister & Lahav 2004| 
software and another red-sequence method (Song et al, m 
prep). This independent cross-checking led to consistent 
redshift estimates in all cases, except for SPT-CL J0615- 
5746, where redshift estimates ranged from about 0.9 to 
1.1 (this cluster was later spectroscopically confirmed at 
z = 0.972); and SPT-CL J0555-6405, where different es- 
timates gave redshifts of 0.27, 0.35, and 0.42. In the 
latter case, because the red sequence appears most pro- 
nounced and u nambiguous at 0.42, we report the redshift 
from the | H10| software only. 

An additional cross-check for two of the clust e rs pre - 
sented here is provided by Menanteau et al.| ( |2010 ), 
who estimated photometric redshifts of 0.75 ± 0.04 for 
SPT-CL J0102-4915 (compare our result of 0.78), and 
0.54 ± 0.05 for SPT-CL J0438-5419 (compare 0.45). 



5.2. Spitzer Space Telescope Imaging 

Spitzer/IRAC imaging is particularly important for the 
confirmation and study of high-redshift SPT clusters, 
such as SPT-CL J2106-5844 at z = 1.132 (Figure [261), 



where the optically faint members are strongly detected 
in the mid-infrared. Three of our catalog clusters were 
observed as part of a larger program to follow up clusters 
identified in the SPT survey. The on-target observations 
consisted of 8 x 100 s and 6 x 30 s dithered exposures at 
3.6 and 4.5 /im, respectively. The deep 3.6 /im observa- 
tions are sensitive to passively evolving cluster galaxies 
down to 0.1 I/* at z = 1.5. The data were reduced ex- 



actly as in Brodwin et al. 



Ashby et at[tlj009| ). B 



after HIO). Photometry wa s calibrated usmg the 
Locus Regression method ("High et al."'2009"). The red- 
sequence redshifts reported in Table [3] were derived using 



2010| ), following the method of 
Ly, we correct for column pull- 
down and residual image effects, mosaic the individual 
exposures, resample to 0^'86 pixels (half the solid angle 
of the native IRAC pixels), and reject cosmic rays. 

5.3. Spectroscopy 

Eleven of the clusters in this work have published spec- 
troscopic redshifts, which we note in Table |3] Using 
the instruments listed in Table [5j we present new spec- 
troscopic redshift measurements on five clusters, four of 
which have no such previously published data. The ro- 
bust biweight location estimator is used to determine the 
cluster spectroscopic redshifts from ensembles of member 
galaxies. 

5.3.1. SPT-CL J0234-5831 and SPT-CL J0254-5857 
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TABLE 3 
Cluster Redshift Data 



Object Name 



Spectroscopy Ref.^ 



# members Zrs Imaging Ref.^ 
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CL 


J0555- 


-6405 
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CL 


J0615- 


-5746 





972 


SPT- 


CL 


J0628- 


-4143 





176 


SPT- 


CL 


J0638- 


-5358 





222 


SPT- 


CL 


J0645- 


-5413 





167 


SPT- 


CL 


J0658- 


-5556 





296 


SPT- 


CL 


J2023- 


-5535 





232 


SPT- 


CL 


J2031- 


-4037 





342 


SPT- 


CL 


J2106- 


-5844 


1 


132 


SPT- 


CL 


J2201- 


-5956 





098 


SPT- 


CL 


J2248- 


-4431 





348 


SPT- 


CL 


J2325- 


-4111 






SPT- 


CL 


J2337- 


-5942 





776 


SPT- 


CL 


J2344- 


-4243 






NOTE. - 


— Spectroscopic and 



raI] ( [T994| 



sp3 

de Grandi et al. 



(1999 



de Grandi et ai.'('iyy9| 
de Grandi et al. (lyW) 
I'ucker et al. ( iy9« ) 
l^onrmger et al. ( :^004b ) 
liolirmger et al. ( 2UU4b ) 
sp4, see 1^'oley et al. (2Ull) 
Struble 6l Kood (1999) 
[j^onrmger et al.[ ( [2UUlt)| > 

sp5; also spl, seelHlOl 



20 



0.40 iml 

0.78 iml 

0.44 im3, im4 

0.53 iml, im4 

0.35 im4 

0.43 im3, im4 

0.52 iml 

0.42 im4 

0.62 iml, im7 

0.45 iml, im4 

0.32 iml 

0.42 iml 

1.0 iml, im6 

0.21 iml 



0.23 im2, im3, im4 

1.17 im5, im6, im7 
im2, im6 

0.37 iml 

0.77 im2, seelmol im7 
0.62 iml 



^ Cross-reference to spectroscopic redshift data from external and internal sources. Internal references 
are denoted "sp#" and refer to Table [s] 

^ Cross-reference to broadband redshift data from external and internal sources. Internal references 
are denoted "im#" and refer to Table ^ 



TABLE 4 
Optical and infrared imagers 



Alias"^ 


Site 


Telescope 


Aperture 
(m) 


Camera 


Filters^ 


Field 


Pixel scale 

n 


iml 


Cerro Tololo 


Blanco 


4 


MOSAIC-II 


griz 


36' X 36' 


0.27 


im2 


Las Campanas 


Magellan/Baade 


6.5 


IMACS f/2 


griz 


27.'4 X 27.'4 


0.200 


im3 


Las Campanas 


Magellan/Clay 


6.5 


LDSS3 


griz 8 


.'3 diam. circle 


0.189 


im4 


Las Campanas 


Swope 


1 


SITe3 


BVRI 


14.'8 X 22.'8 


0.435 


im5 


Par anal 


VLT 


8.2 


FORS2 


I 


6.'8 X 6.'8 


0.25 


im6 


Cerro Tololo 


Blanco 


4 


NEWFIRM 


JKs 


28' X 28' 


0.4 


im7 




Spitzer Space Telescope 


0.85 


IRAC 


[3.6] [4.5] 


5.'2 X 5.'2 


1.2 


Note. — The optical and infrared cameras used. The choice of facilities and filters for any given cluster was typically optimized 
according to our best redshift estimate prior to observation. Not all imagers, nor all the listed filters, were used on each cluster. 
^ Shorthand alias used in Table [s] 

^ The filters we used, which were in general a subset of all of those available. We did not typically use all listed filters on each 



cluster. 



TABLE 5 

Optical and infrared spectrographs 



Alias^ 


Site 


Telescope 


Aperture 
(m) 


Camera 


Mode 


spl 


Las Campanas 


Magellan /Clay 


6.5 


LDSS3 


longslit 


sp2 


Las Campanas 


Magellan/Baade 


6.5 


IMACS 


GISMO 


sp3 


Las Campanas 


Magellan/Baade 


6.5 


IMACS 


longslit 


sp4 


Par anal 


VLT 


8.2 


FORS2 


MOS 


sp5 


Cerro Pachon 


Gemini South 


8.1 


CMOS 


MOS 



Note. — The spectrographs used. 
^ Shorthand alias used in Table [s] 
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In a procedure similar to Brodwin et al. (2010), multi- 
slit spectroscopic observations were acquired on the 6.5- 
meter Baade Magellan telescope on UT 2010 October 8 
for SPT-CL J0234-5831 and SPT-CL J0254-5857. Mea- 
surements were made using the Gladders Image- Slicing 
Multislit Option (GISMO, Gladders et al. in prep.) mod- 
ule on the Magehan/IMACS spectrograph. GISMO op- 
tically remaps the central region of the IMACS field of 
view (roughly 3.5' x 3.2') to sixteen evenly-spaced regions 
of the focal plane, allowing for a large density of slit lets 
in the cluster core while minimizing trace overlaps on the 
CCD. 

In designing the multislit mask, galaxies were assigned 
a weight proportional to their r-band brightness and ad- 
justed for their position in color space with respect to 
a manually-selected red sequence. The f/4 camera, the 
300 1/mm grating and the zl430 — 675 filter were used. 
Each cluster was observed with three 30-min exposures 
of one mask in good seeing (~ 0.6"). 

The COSMOS reduction package was used for stan- 
dard CCD processing, resulting in wavelength-calibrated 
2D spectra. The ID spectra were then extracted from 
the sum of the reduced data. Secure redshifts were ob- 
tained for 21 member galaxies of SPT-CL J0234-5831, 
and 32 member galaxies of SPT-CL J0254-5857. 

5.3.2. SPT-CL J0615-5746 

Longslit spectroscopy of SPT-CL J0615-5746 was per- 
formed on UT 2011 March 8, also with the IMACS spec- 
trograph on the Baade Magellan telescope. The longslit 
was aligned across several objects and yielded clear red- 
shifts for the BCG and a second cluster member. The 
reported redshift is that of the BCG. 

5.3.3. SPT-CL J2106-58U 



We refer the reader to Foley et al. (2011) for a de- 



tailed discription of spectroscopic measurements of SPT- 
CL J2106-5844. In short, the redshift given in Tablejs] 
is derived from 18 member galaxies using VLT/F0RS2 
and Magellan/IMACS-GISMO. 

5.3.4. SPT-CL J2337-5942 
The redshift of SPT-CL J2337-5942 reported in Ta- 



ble [6| 



^spec 



0.776, is from combined measurements of 



19 cluster members using CMOS on the 8.1 m Gem- 
ini South telescope and 2 members using the Magel- 
lan/LDSS3 longslit — one of which overlaps with a CMOS 
member — for a total of 20 cluster members. The LDSS3 
data were described in detail in HIO^ where ^spec = 0.781 
was reported from the two members. 

For the new CMOS observations we are presenting 
here, galaxies with r — i color consistent with a cluster 
red- sequence at z = 0.77, and having non-stellar PSFs 
in the Gemini i-band pre-image, were used to populate 
two masks. A total of 31 galaxies were observed for three 
hours with the R150_G5326 grism and the GG455-G0329 
filter. The IRAF Gemini reduction package was used for 
standard CCD processing. The wavelength-calibrated 
2D spectra were sky-subtracted using an in-house rou- 
tine, after which the ID spectra were extracted from the 
coadded 2D spectra. Secure redshifts for the 19 cluster 
members were obtained from the CMOS masks. All the 
spectra are from early- type galaxies, often exhibiting a 
very strong Ca H&K absorption feature. 



6. X-RAY DATA 

For each cluster in Table [6) we searched the ROSAT 
data archive for possible X-ra y counterparts, in clud- 
ing the REFLEX catalog (Bo hringer et al.||2004a1), the 
ROSAT A ll-Sky Survey Bright Source Catalog (|Voges 



et al.|1999[ ), and the ROSAT All-Sky Survey Faint Source 
Catalog^ We also used data from the ROSAT Ah 
Sky Survey and pointed Position Sensitive Proportional 
Counter (PSPC) observations to measure the X-ray fiux 
for the SPT clusters in the "hard" 0.6-2 keV energy band, 
since the signal to noise is generally bet ter in this band 
than in the full ROSAT energy range (Vikhlinin et al. 
|1998| ). To determine the source counts, we used a source 
radius corresponding to 2 Mpc, excluded any sources not 
associated with the cluster emission and used a nearby 
region for measuring the X-ray background. Significant 
detections were found for 25 of the 26 clusters in the 
ROSAT observations. The highest-redshift cluster, SPT- 
CL J2106-5844, was not detected by ROSAT, but has 
been detected in pointed observations with Chandra ( jFo^ 
ley et al. 2011). 

To determine the cluster fiux and luminosity, we used 
PIMMS^^ to determine the cluster unabsorbed fiux in the 
observer's frame from the ROSAT "hard" band observa- 
tions. We then used the XSPEC^^ "fiux" and "lumin" 
functions to determine the cluster luminosity in the clus- 
ter rest frame, which we report in Table [6j We report 
each in the 0.5-2.0 keV band because of this band's rela- 
tive insensitivity to the assumed X-ray temperature. For 
example, the inferred fiux changes by ^2% when assum- 
ing a range of gas temperatures between 6 to 10 keV. For 
two clusters, SPT-CL J2106-5944 and SPT-CL J2337- 
5942, we give the fiu x and luminos i ty me asure d by Chan- 
dra and rep orted in Foley et al. (2011) and Andersson 
et al. (|2010"), respectively. 



TiTFigure [4j we plot the X-ray luminosity, Lx, and 
the SPT measured massfconverted from M2oo(Pmean) to 
M5oo(Pcrit)) from Table Im We plot only statistical uncer- 
tainties for both luminosity and mass. However, we note 
that we have ignored several important systematic uncer- 
tainties in the ROSAT X-ray luminosity measurements, 
including the effects of unresolved point sources and cool- 
ing cores. These can add biases and additional scatter to 
the X-ray measurements, however cannot observationally 
be accounted for because of the relatively large ROSAT 
beamsize. Similar phenomena also affect the SZ measure- 
ments, however they are accounted for statistically in the 
significance estimate and the scatter in the significance- 
mass relation, as described in Section |4.1[ The latter 
dominates the statistical uncertainty of the SZ mass esti- 
mates. In Figure [4j we have assumed the best- fit redshift 
evolution in t he Lx-M relation measured by Vikhlinin] 
where Lx cx ME{zy-^^. This is slightly 



et al. 



( |2009aD , 

ditterent than the self-similar expectation that predicts 
an evolution of E{z)'^ for the luminosity in the 0.5-2.0 
keV band and an evolution of ^ E{zy^^ for the bolomet- 
ric luminosity, which has an additional weak dependence 
of luminosity with temperature. 
In Figure [4j we also show the best-fit Lx-M rela- 

"^■^ http://heasarc.gsfc. nasa.gov/W3Browse/rQsat/rass fsc. html| 
http://cxc.harvard.edu/toolkit/pimms.jsp 
http://heasarc.nasa.gov/xanadu/xspec/xspecll/index.html 
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SZ (10^^ M3„|/hJ 



Fig. 4. — The X-ray luminosity and SZ inferred masses 
-^50o(Pcrit) for our cluster sample. We plot statistical uncertainties 
only, and note that the statistical uncertainty of the SZ mass esti- 
mate is limited by the assumed scatter in the SZ significance-mass 
relation. Clusters from the shallow fields are in blue, and clusters 
from the deep fields are in red. We also show the best-fit relations 
of Pratt et al. (2009) (dotted), Vikhlinin et al. (2009a) (dash-dot), 
and Mantz et al., ( |2010, ) (dashed). 



tions from Pratt et al. (2009), Vikhlinin et al. (2009a 



and Mantz et al. ( |20T0 ), to co mpare our results with 
other X-ray cluster studies. BothlPratt et al. (2009) and 
Mantz et al. (2010) assume evolution consistent with the 
selt-similar relation for bolometric luminosity. For each 
result, we use their published normalization and slope, 
but have assumed the redshift evolution measured by 
Vikhlinin et al. (2009a). Relative to self-similar evolu- 



tion, this would cause a < 7% change in the normaliza- 
tion for the typical redshif t range of their cluster samples 



{z < 0.3). [Mantz et al 
nosities in the rest frame 



(|2010| also quote their lumi- 
lA keV band, which we 
have converted to 0.5-2.0 keV by dividing by a factor 
of 1.61. This factor is appropriate for a cluster with a 
8 keV electron temperature, and varies negligibly with 
temperature. We expect that both approximations will 
not significantly change the normalization of either re- 
lation. As can be seen in Figure |4j there is a fairly 
large spread in the published Lx-M relations, although 
the methods and samples differ between each work. For 
example, they each use cluster samples with somewhat 
different mass and redshift ranges, and account for the 
effects of Eddington bias differently. Regardless, we con- 
sider the agreement with other cluster samples reason- 
able, and also significant given the unique SZ selection 
of the clusters in this work. The X-ray selected cluster 
samples are of generally lower redshift, less massive, and 
have been corrected for Malmquist bias from the X-ray 
fiux selection, a bias that is completely absent for SZ se- 
lected samples. This work confirms that the SPT cluster 
sample consists of very massive clusters, which qualita- 
tively follow the Lx-M relation measured from other 
X-ray selected cluster samples. 

7. DISCUSSION 



The 26 highest-significance SZ-selected clusters from 
the 2500 deg^ SPT survey (see Table ^ include all the 
most massive galaxy clusters in this region of the sky, 
independent of the cluster redshift. These exceedingly 
rare systems populate the high end of the mass function 
at each redshift, and predictions for the characteristics of 
this population are sensitive to the details of the assumed 
cosmological model. An interesting first step in using the 
clusters presented in this work to constrain cosmology is 
to ask whether their distribution in mass and redshift is 
consistent with the predictions of the standard ACDM 
cosmological model. We investigate this quest ion two 
ways. First, we use the framework of Mortonson et al. 



(2011, hereafter Mil) and the fitting functions they pro- 
vide to ask whether the existence of any single cluster in 
our sample is in significant tension with ACDM. We then 
fit all available cosmological data including this new clus- 
ter sample to two different cosmological models — namely, 
standard ACDM and a single-parameter extension allow- 
ing for non- Gaussian initial conditions — and see if the 
data prefer the non-standard model. 

7.1. Single- duster Tests 

Mil have published fitting functions that allow us to 
answer the question: Is this one cluster in significant 
tension with ACDM? In Figure [5j we plot the mass vs. 
redshift for all 26 clusters and overplot exclusion curves 
from Mil. As explained in Appendix C of Mil, the 
mass for a given cluster that is appropriate to compare 
to their exclusion curves is not precisely the best pos- 
terior estimate of that cluster's mass. The masses plot- 
ted in Figure |5] are calculated using Equation C3 from 
Mil, using the conditional probability P(^|M) to esti- 
ma te their Mphs and a ]n m and using the local slope of 
the [Tin ker et al I (]2008|) mass function around the value 
of Mobs foi" 7- Fhe error bars on the masses in Figure [5] 
are calculated by setting the fractional error (i.e., the er- 
ror in In M) equal to the fractional error of the posterior 
mass estimates in Table [U We have confirmed that this 
is an excellent approximation to the full probability dis- 
tribution of the Mil-appropriate masses for the cluster 
sample in this work. 

The two highest exclusion curves overplotted in Figure 
[5] represent the mass and redshift above which an indi- 
vidual cluster would be less than 5% likely to be found 
in a given survey region in 95% of the ACDM parame- 
ter probability distribution. We plot one exclusion curve 
for the least likely cluster allowed in a 2500 deg^ survey 
and one curve for the least likely cluster allowed in the 
entire sky. It is clear that, according to the formalism of 
Mil, no cluster in our sample is individually in strong 
tension with ACDM — a conclusion that would still hold 
if we appli ed t he naive scaling factor of 1/0.78 discussed 
in Section [m to all the cluster masses. 



This result can be compared to the result of [Foley et al. 



([20TT|, in which the single cluster SPT-CL J2106-5844 
is found to be less than 5% likely to exist in the 2500 
deg^ SPT survey region in 32% of the ACDM parameter 
probability distribution. There are some differences in 
the two ana l yses, the most important of which is that 



Foley et al. (2011) use a mass estimate that combines 
and X-ray data, whereas this work onl y reports an 
SZ- der ived mass. The central value of the 'Foley et al.[ 
(2011) combined SZ/X-ray mass estimate is 30% higher 
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Fig. 5. — An Mil-style plot showing the mass M2oo(pmean) and 
redshift of the clusters presented in this paper. The masses in 
this plot are slightly (~5%) higher than those presented in Table 
[g] due to the slightly different treatment of mass bias appropriate 
to the Mil calculation (see text and Mil for details). Some of 
the most extreme objects in the catalog are annotated with the 
R.A. portion of their object name. The red solid line shows the 
mass above which a cluster at a given redshift is less than 5% 
likely to be found in the 2500 deg^ SPT survey region in 95% 
of the ACDM parameter probability distribution. The black dot- 
dashed line shows the analogous limit for the full sky. The blue 
open data point (redshift slightly offset for clarity) denotes the 
mass estimate for SPT-CL J2106-5844 from combined X-ray and 
SZ measurements in Foley et al. (2011). That work concludes that 
this cluster is less than 57o likely in 32% of the ACDM parameter 
probability distribution, and we show the corresponding Mil p = 
32% limiting mass vs. redshift as the dashed blue line. 



than the central value of t he SZ- derive d mass reported 
here. We have included the [Foley et al. (2011) combined 
mass as a point in Figure [Sj and we have also plotted the 
Mil exclusion curve corresponding to < 5% likelihood of 
finding a cluster in the SPT survey in 32% of parameter 
probability. As expected from the result in Foley et al. 
f2011), the p = 32% exclusion curve nearly intersects the 
central SPT-CL J2106-5844 mass value from that work 
(adjusted appropriately for the Mil plot). 



7.2. Extensions to ACDM 

While no individual cluster lies above either p = 95% 
exclusion line in Figure [5j there are several which come 
reasonably close. One might imagine that the collective 
"unlikelihood" of these clusters could indicate the need 
to go beyond the standard ACDM cosmological model. 
The most straightforward extension to ACDM that could 
explain an excess of massive clusters (including ones at 
high redshift) is the possibility of a non-Gaussian compo- 
nent to the primordial density perturbations. Different 
models of infiation predict different levels and types of 
non-Gaussianity (e.g., Bartolo et aL||2QQ4 ), with the size 
of the leading-order non-Gaussian term described by the 
parameter /nl- We have included the mass and redshift 
distribution of the clusters presented here in a cosmolog- 
ical likelihood calculation with and without /nl as a free 
parameter. The likelihood calculation was implemented 



as in VIO, with the effect of /nl o n cluster abundanc e 
added following the prescription of Dalai et al. ( |2QQ8 ). 
To simplify the selection function and mass scaling part 
of the calculation, the preview-depth relation (see Table 

?was used for all clusters, and preview-depth values of 
were estimated for the full-depth clusters by making 
coadded maps of only one ninth of the observations and 
running the cluster finder on these maps. As in VIO, the 
scaling relation is not expected to capture the correct be- 
havior at low redshift; as in VIO, we exclude this regime 
by applying a hard cut z > 0.3 in this analysis. 

The preferred value of /nl in the extended model is 
consistent with zero {/nl = 30 ± 450 at 68% confi- 
dence). We note that this is a significantly weake r con- 
straint than that found by [Komatsu et al. (2011) using 
the CMB bispectrum as measured in the WMAFy data, 
but that the two results are consistent with each other 
and with f]^r. = 0. This is in tension with the recent 



results of [Cayo n et al.) ( |2010D , [Hoyle et alT] ( |2010D , and 
Enqvist et al.| ( |2uTol ), who found signiticant evidence for 
non-zero /atl based on other high redshift galaxy clus- 
ters using a different statistical technique. In contrast 
to those works, this analysis uses a likelihood analysis 
over the full range of mass and redshift space including 
the SPT selection function, marginalizing over scaling 
relation uncertainties. This approach naturally incor- 
porates information about both the mass and redshift 
distribution and the total number of clusters. Our /nl 
constraints do not change appreciably and are still fully 
consistent with /tvl = if we re-run the analysis using 
a prior on the ^-M relation that incorporates the scaling 



between SZ and X-ray masses discussed in Section 4.1 



8. CONCLUSIONS 

We have presented a sample of the most massive galaxy 
clusters in a 2500 deg^ region of the sky, selected via their 
SZ signature in SPT observations. These 26 clusters are 
selected from hundreds of SPT cluster candidates on the 
basis of their SZ detection significance, which has been 
shown in previous SPT analyses to correlate tightly with 
cluster mass (VIO, Anderss on et al.||2010[ ). As expected 
from their high SZ significance, each one of these objects 
shows a strong overdensity of similarly colored galaxies in 
optical and/or infrared data, and the X-ray luminosity of 
these systems (as estimated from archival and newly col- 
lected data) is consistent with their SZ-derived masses. 

We measure (or collect from the literature) 
photometric — and, in some cases, spectroscopic — 
redshifts for these 26 clusters. The cluster sample 
includes several newly discovered high-redshift systems, 
significantly increasing the total number of known 
galaxy clusters with masses of M2oo(Pmean) ^ 10^^ Mq 
and redshifts z > 0.5. In addition to being interesting 
targets for studies of cluster physics and galaxy for- 
mation in the densest environments in the Universe, 
these massive, high-redshift clusters allow us to test 
the standard ACDM cosmological model with Gaussian 
initial conditions. 

We test whether the most extreme (in mass or red- 
shift) individual clusters pose a challenge to ACDM by 
applying the formalism presented by Mil. No single clus- 
ter is in significant tension with the ACDM model, with 
the caveat that current SZ mass estimates are uncer- 
tain at the 30% level. Improved mass estimates would 
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strengthen any test of the ACDM model. We also exam- 
ine constraints on primordial non-Gaussianity from the 
cluster sample. The data show no preference for non- 
Gaussianity as parametrized hj /nl- At 68% confidence, 
we find fNL = 20 =b 450. 

The sample of clusters found with the SPT, the most 
massive of which are presented in this work, is comple- 
mentary to the sample expected to be found with the 
Planck satellite. Planck will find very massive galaxy 
clusters over the entire sky — 16 times more area than 
the SPT survey. However, clusters above z ^ 0.5 will 
have a typical angular size of 1', meaning that Planck^ s 
sensitivity to clusters will fall off at high redshift due to 
its larger beam (5' FWHM at 150 GHz vs. I' for SPT). 
The clusters found by Planck will thus be distributed 
towards significantly lower redshift than those found by 
the SPT. 

The SPT observations of the 2500 deg^ of sky used in 
this work will be completed to the final survey depth in 
the Austral winter of 2011. Extrapolating from current 
survey yields, the complete SPT SZ survey will contain 
roughly 750 galaxy cluster candidates at a detection sig- 
nificance of ^ > 4.5 (and over 400 at ^ > 5), the vast 
majority of which will correspond to real, massive clus- 
ters. This unique, nearly mass-limited cluster sample will 
offer an unprecedented opportunity to test the ACDM 
cosmological model and the properties of dark energy. 

Facilities: Blanco (NEWFIRM), Blanco (MOSAIC), 
CXO (ACTS), Gemini-S (CMOS), Magellan: Baade 
(IMACS), MagellaniClay (LDSS3), Spitzer (IRAC), 
South Pole Telescope, VLTiAntu (F0RS2) 
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APPENDIX 



SZ AND OPTICAL/INFRARED IMAGES 



Figures 6-31 show SZ detection significance maps {left panels) and optical and infrared images {right panels) of the 
clusters. In all images, north is up, east is left. The SZ-only insets subtend 12 arcminutes on a side. The mapping 
between color and SZ significance ^ is different in all SZ thumbnails, spanning the full range of SZ pixel values in the 
region of sky shown. The peak value in each thumbnail is equal to the quoted SZ detection significance in Table [6) 
Contours denote significance values of (—8, —4, —2, 0, 2, 4, 8, 16, 32) in all thumbnails. Contours are dashed where ^ is 
negative, and solid where ^ is positive. The negative lobes around some of the most significantly detected clusters in 
the SZ images are due to the filtering of the time-ordered data and the maps 

The optical/infrared images have the same contours as their corresponding SZ thumbnail overlaid. False-color 
composites are presented for clusters where multiband imaging is available, either from our own observations or from 
public archives. Otherwise, black-and-white images are shown. 




Fig. 7. SPT-CL J0102-4915, also known as ACT-CL J0102-4915, at Zrs = 0.78. Blanco/MOSAIC-II irg images are shown in the 
optical/infrared panel. 




Fig. 9. — SPT-CL J0234-5831 at Zspec = 0.415. Magellan/LDSS3 zrg images are shown in the optical/infrared panel. 



41.05 41.00 40.95 40.90 40.85 40.80 
R.A. 



40.98 40.96 40.94 40.92 40.90 40.88 40.86 40.84 

R.A. 



Fig. 10. SPT-CL J0243-4833 at Zrs = 0.44. Blanco/MOSAIC-II irg images are shown in the optical/ infrared panel. 





Fig. 13. — SPT-CL J0304-4401 at Zrs = 0.52. Blanco/MOSAIC-II irg images are shown in the optical/ infrared panel. 




Fig. 15. SPT-CL J0417-4748 at ^rs = 0.62. Spitzer/IRAC [3.6] and Blanco/MOSAIC-II ig images are shown in the optical/ infrared 
panel. 




Fig. 17. SPT-CL J0549-6204 at Zrs = 0.32. Blanco/MOSAIC-II irg images are shown in the optical/ infrared panel. 




Fig. 19. SPT-CL J0615-5746 at Zspec = 0.972. Blanco/NEWFIRM Ks and Blanco/MOSAIC-II rg images are shown in the opti- 
cal/infrared panel. 




Fig. 21. SPT-CL J0638-5358, also known as ABELL S0592, RXCJ0638. 7-5358, and ACT-CL J0638-5358, at Zspec = 0.222. 
Spitzer/IRAC [4.5] [3.6] and Gemini/GMOS r images are shown in the optical/infrared panel. 




Fig. 23.— SPT-CL J0658-5556, also known as the Bullet Cluster, RXCJ0658. 5-5556, and ACT-CL J0658-5557, at Zspec = 0.296. MPG- 
ESO/WFI R and V images are shown in the optical/infrared panel. 




Fig. 25.— SPT-CL J2031-4037, also known as RXC J2031.8-4037, at Zsp, 
optical/infrared panel. 



= 0.342. Magellan/LDSS3 irg images are shown in the 




Fig. 27.— SPT-CL J2201-5956, also known as Abell 3827 and RXCJ2201. 9-5956, at Zspec = 0.098. IMACS f/2 irg images are shown in 
the optical/infrared panel. This detection is at the eastern edge of the survey field. 




Fig. 29.— SPT-CL J2325-4111, also known as Abell S1121, at Zrs = 0.37. Blanco/MOSAIC-II irg images are shown in the optical/infrared 
panel. 




Fig. 31. SPT-CL J2344-4243 at Zrs = 0.62. Blanco/MOSAIC-II irg images are shown in the optical/infrared panel. 



